TRAFFIC DATA AND THEIR IMPLICATIONS FOR 
CONSISTENT TRAFFIC FLOW MODELING 



Dirk Helbing*'^ 

* //. Institute of Theoretical Physics, University of Stuttgart, 
Pfaffenwaldring 57/111, 70550 Stuttgart, Germany, 
www: http://www.uni-stuttgart.de/UNIuser/thphys/helhing.html 



Abstract: The paper analyzes traffic data with respect to certain aspects which are 
relevant for traffic flow modeling as well as the calibration of model parameters and 
functions. Apart from the dynamic velocity distribution, the density-dependence and 
the temporal evolution of various, partly lane-specific quantities is investigated. The 
results are well compatible with recent macroscopic traffic ffow models which have 
been derived from the dynamics of driver-vehicle units. These have also solved the 
inconsistencies, which previous models have been criticized for. 
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1. INTRODUCTION 

Efficient traffic flow is very important for a 
high standard of living. However, during the last 
decades the traffic volume has continuously in- 
creased, so that the capacity of the road system 
is almost reached in many countries, especially 
in conurbations. Therefore, methods for traffic 
flow optimization are called for. To develop these, 
numerous traffic simulation models have been pro- 
posed. 

This paper will focus on macroscopic (fluid- 
dynamic) models for the vehicle density and the 
average velocity. It addresses essentially two ques- 
tions: 1. Which phenomena and functional rela- 
tions should be reproduced by the models? Re- 
lated results are presented in the next section. 
2. Which one of the existing traffic ffow models 
should be used to simulate traffic dynamics? This 
is the topic of Sect. 3. In particular, it is discussed 
which consistency requirements should be met by 
realistic models. 



^ The author is grateful to Henk Taale and the Ministry 
of Transport, Public Works and Water Management for 
supplying the freeway data. 



2. DATA EVALUATION 

During the last year single vehicle data of the 
Dutch two-lane Autobahn A9 from Haarlem to 
Amsterdam have been evaluated with respect to 
various criteria which are relevant for traffic flow 
modeling. The quantities measured by induction 
loops at discrete places x along the freeway in- 
clude the passage times ta{x), velocities Va{x), 
vehicle lengths la{x), and lanes ia{x) of the single 
vehicles a. These allow the determination of var- 
ious aggregate (so-called macroscopic) quantities. 

The results of the standard evaluation methods 
have already been presented, compared and dis- 
cussed elsewhere (Helbing, 1997). They base on 
averages over a fixed time interval T, giving large 
statistical variations of the velocity moments 

^ ' t~-T/2<ta<t+T/2 

where the traffic flow Q{x, t) = N{x, t)/T is small, 
since the number N{x,t) of vehicles passing the 
intersection at place x between times t — T/2 and 
t + T/2 is smaU, then. 



The following considerations will focus on an eval- 
uation method with a sampling error that is inde- 
pendent of the traffic volume. For this purpose, 
one averages over a fixed number N = 100 of 
vehicles, implying a variation of the time interval 
T{x,t) of data collection. Then, the definition 
of the space- and time-dependent traffic fiow is 
Q{x, t) = N/T{x, t) with 



t = 



and the velocity moments are given by 
1 
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ao<a<ao+N 



This allows the calculation of the average ve- 
locity V{x,t) = {v), the density p{x,t) = 
Q{x,t)/V{x,t), the velocity variance 

e{x,t) = {[v-Vix,t)r) = {v')-{vf, (4) 
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The definitions for the single lanes i are analogous 
to the above ones for the total cross section. 

Whereas the variance is a measure for the 
breadth of the velocity distribution P{v; x, t), the 
skewness is a measure for its asymmetry, and 
the kurtosis for its flatness. Most previous studies 
have restricted to the evaluation of the velocity- 
density relation Ve{p) (Fig. 1) and the fundamen- 
tal diagram Qe(p) (Fig- 2). 

However, for some theoretical considerations, it is 
essential to know the form of the velocity distribu- 
tion in dependence of density. In agreement with 
previous studies it was found that a normal dis- 
tribution is a good approximation for the average 
velocity distribution (Fig. 3). Recent numerical 
solutions for the stationary velocity distribution 
of an improved kinetic traffic model yield similar 
results (Fig. 4). 

Due to the limited amount of data, the time- 
dependent velocity distribution is hard to obtain. 
Therefore, it has been investigated whether the 
time-dependent skewness (Fig. 5) and the kurto- 
sis (Fig. 6) vanish in agreement with a normal 
distribution. The deviation of the skewness from 
zero is indeed neither systematic nor significant. 
The kurtosis is also very small but, on average, 
somewhat smaller than zero. Since this systematic 
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Fig. 1. Empirical velocity- density data (•) for the 
Dutch two-lane freeway A9 with a speed limit 
of 120km/h. The fit function Ve{p) {—) at 
high densities is reconstructed by means of 
theoretical traffic flow relations (cf. Sect. 3). 
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Fig. 2. Empirical flow-density data (•) and fit of 
the fundamental diagram Qe(p) ( — )• 
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Fig. 3. Comparison of empirical velocity distri- 
butions at different densities ( — ) with fre- 
quency polygons of grouped normal distri- 
butions having the same mean value and 
variance ( — ). A significant deviation of the 
empirical relations from the respective dis- 
crete normal approximations is only found at 
p = 40 vehicles/km lane, where the averaging 
interval T may have been too long due to 
rapid stop-and-go waves (cf. the mysterious 
"knee" at p « 40 vehicles/km lane in Fig. 7). 
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Fig. 4. Stationary velocity distributions at dif- 
ferent vehicle densities obtained from a im- 
proved kinetic traffic model (reproduction by 
kind permission of C. Wagner). 

deviation from zero can be shown to vanish on the 
left-hand lane, it is probably an effect of the finite 
truck fraction driving slower than other vehicles. 

In summary, the velocities are almost normally 
distributed, even in the dynamic case. Neverthe- 
less, if one averages over too large time intervals 
T, bimodal distributions may result in cases of 
strong variations of traffic flow. 
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Fig. 5. The empirical skewness (•) is approxi- 
mately zero most of the time. 

Moreover, it is interesting to know the variance- 
density relation (Fig. 7). Looking at the temporal 
evolution of the variance, one detects significant 
peaks when the average velocity breaks down or 
recovers (Fig. 8). Dctaillcd mathematical consid- 
erations show (Helbing, 1997), that it can be re- 
constructed from the temporal course of the vehi- 
cle density p{x, t) via the variance-density relation 
0e(/3), if one takes into account a positive correc- 
tion term which is caused by averaging over finite 
time periods T. The theoretical relation reads 



4 



dV 



(7) 



Investigating the temporal fluctuations of the av- 
erage velocity V{x,t), one finds a fiat power spec- 



Fig. 6. The empirical kurtosis shows a small sys- 
tematic deviation from zero due the finite 
truck fraction. Nevertheless, it is almost zero. 
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Fig. 7. Standard deviation •y/8(x, t) of vehicle 
velocities in dependence of density p{x. t) (■) 
and corresponding fit function y/OJJj) ( — ). 
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Fig. 8. The standard deviation ^yQ{x,t) of vehicle 
velocities ( — ) has maxima, when the average 
velocity V{x, t) ( — ) changes considerably. It 
cannot be described by the variance-density 
relation ^^/Qf,[p{x,t)) alone ( ), but to- 
gether with an additional positive term which 
arises from the evaluation procedure (•••). 



trum Vi^u), i.e. almost the same intensity of the 
different fluctuation frequencies, corresponding to 
a so-called "white noise" (Fig. 9). However, at 
small frequencies v (large oscillation periods l/i/) 
the power law 

V{u) = Cv-^ (8) 

appears. Simular results are obtained for the fluc- 
tuations of the density p{x,t) or the variance 
Q{x,t). The power law behavior is most distinct 
near to entrances, but only weakly visible at 
undisturbed sections of the freeway. 
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Fig. 9. Power spectrum of the average velocity's 
temporal course ( : logC — 5\ogv, 5 € 

{0,2}). 

It has also been found that the dynamics on 
neighboring lanes is strongly correlated (Fig. 10). 
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Fig. 10. Comparison of the temporal course of the 
vehicle density Pi{x,t) on the left lane ( — ) 

with the right lane ( ) and of the average 

velocity Vi{x, t) on the left lane ( — ) with the 
right lane (•••)• 

This justifies the common practice to model the 
dynamics of the total cross section of the road in 
an overall manner, treating it like a single lane 
with higher capacity and possibilities of overtak- 
ing (instead of explicitly describing the dynamics 
of all the lanes, including their interactions by 
overtaking and lane-changing maneuvers). Above 



an average density of 35 vehicles per kilometer and 
lane the velocity on the left lane is the same as on 
the right lane, but the left lane is more crowded 
than the right one. The variance on neighboring 
lanes behaves almost identical (Helbing, 1997). 

The maximum flow is 2200 vehicles per hour and 
lane, and it is reached at a density of 30 vehicles 
per kilometer and lane. Above this density, traffic 
flow breaks down and stop-and-go traffic develops. 
However, the hysteresis curve in Fig. 11 indicates 
that traffic flow is already unstable at densities of 
about 12 vehicles per kilometer and lane. 
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Fig. 11. Comparison of the fundamental diagram 

Qe{p) ( , cf. Fig. 2) with the temporal 

evolution of traffic flow Q{x,t) ( — ). 

Although the average velocity V{x,t) recovers 
after a bottleneck (at x = km), the amplitude 
of the emerging stop-and-go waves becomes larger 
and larger (Fig. 12). At the same time their wave 
length grows, indicating a merging of traffic jams. 
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Fig. 12. Temporal evolution of the average veloc- 
ity V{x,t) at successive cross section of the 

freeway {■ ■ ■: x = km; : x = 0.95 km; — : 

X = 2.15km; x = 4.15km). 

3. CONSISTENT TRAFFIC MODELING 

For a realistic description of traffic dynamics, it is 
sufficient to restrict the model to the continuity 
equation 



dp ^ d{pV) ^ ^ 
dt dx 



(9) 



for the vehicle density p{x, t) and a suitable dy- 
namic equation for the average velocity V{x,t). 
The velocity equation is necessary for the de- 
lineation of emergent traffic jams and stop-and- 
go traffic. The dynamics of the variance can be 
reconstructed from the temporal evolution of the 
density and average velocity according to Eq. (7). 
Due to the strong correlation between neighboring 
lanes, an overall description of the total road cross 
section is possible as long as the number of lanes 
does not change. Nevertheless, models for the 
interaction between neighboring lanes have also 
been built and simulated (Hclbing, 1997). 

It can be shown that the velocity equations of 
most previous macroscopic traffic models (in their 
continuous version) can be written in the form 

dV dV IdV 1 

The main difference is the specification of the 

traffic pressure V, the relaxation time r and the 
velocity-density relation Ve{p). For example, the 
Lighthill-Whitham model results in the limit r —> 
0. Payne's and Papageorgiou's model is obtained 
for V = -Ve/2T. For dV/dp = -p/[2t{p + 
K)]dVe/dp one ends up at Cremer's model. In 
the model of Phillips there is V = pQ, where 
& denotes the velocity variance. The model of 
Kiihne, Kerner and Konhauser (Fig. 13) results 
for V = p<do — rjdV/dx, where 9o is a positive 
constant and rj a viscosity coefficient. 
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Fig. 13. The representative streamlines related 
to the simulation of the Kerner-Konhauser 
model indicate the average velocity by their 
slope and the vehicle density by the line den- 
sity. Starting with an almost homogeneous 
traffic situation with a density of 40 vehi- 
cles per kilometer and lane, density clusters 
(which correspond to traffic jams) develop in 
the course of time. These merge into greater 
clusters in agreement with Fig. 12. 



All these models have proved their value in many 
applications, and it is hard to decide which is 
the best one. Nevertheless, recently these models 
have seriously been criticized by Daganzo and 
others. Therefore, a new macroscopic model has 
been derived from Boltzmann-like kinetic traffic 
equations, basing on a microscopic model of the 
dynamics of driver- vehicle units. The kinetic traf- 
fic model is related to Paveri-Fontana's one, but 
it additionally takes into account vehicular space 
requirements (by applying the Enskog-formalism 
for dense gases) and velocity fluctuations due to 
imperfect driving (Helbing, 1997). 

The calculation of the corresponding macroscopic 
traffic equations is very laborious. For the equilib- 
rium velocity one finds the relation 



V%p,V,&) = Vo~j 



Tip)[i - pip)]pe 



P/Pr, 



pTrV 



,(11) 



where Vq is the average desired velocity, p{p) the 
probability of immediate overtaking, Pmax ~ 

160 

vehicles per kilometer and lane the maximum ve- 
hicle density, and Tr ~ 0.8 s the reaction time. The 
equilibrium variance 0'^(p, V, 6) = A{p){y'^ + 6), 
where A{p) denotes the strength of velocity fluc- 
tuations, vanishes with the average velocity 
as required. Utilizing these relations, it is pos- 
sible to reconstruct the velocity-density relation 
Ve{p) = F^(/9,y^e'=) (Fig. 1) and the variance- 
density relation Qe{p) = <ci<'{p,V<' (Fig. 7) 
at high densities p (Helbing, 1997). Applying the 
equilibrium approximation 



e(a;, t) « e''{p{x, t),V{x, t), e%x, t)) 
_ A{p{x,t))[V{x,t)]^ 
1-A{p{x,t)) ' 



(12) 



one finally arrives at the corrected velocity equa- 
tion 



dV_ 

at 



V 



dV 
dx 



IdVdp ^dV ^d'^p 
p dp dx dx dx"^ 



n d'^V 1 

Comparing this with (10), one obviously obtains 
the additional terms adV/dx and —hd'^p/d:r?. 
Moreover, V, rj, V''\ a, and b are not only functions 
of the density p anymore, but they also depend on 
the average velocity V. The corresponding func- 
tional relationships can be analytically derived, 
and the model parameters have been determined 
from empirical data. 

It can be shown that the model automatically 
meets all consistency criteria: 1. The equilibrium 
velocity V'^{p,V) decreases monotonically with 
growing density p and vanishes at the maximum 
vehicle density Pmax- 2. Since dV/dp and b vanish 



for V = 0, the average velocity V cannot be- 
come negative. 3. Since dV/dp is non-negative 
(Fig. 14), the traffic pressure monotonically in- 
creases with the density p. As a consequence, 
vehicles will not accelerate in direction of grow- 
ing density {dp/dx > 0). 4. The traffic pressure 
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Fig. 14. Representation of the density-gradients 
of the idealized traffic pressure V = p&eip) 
of point-like vehicles ( — ) and the corrected 
pressure relation ( — ) which takes into ac- 
count their finite space requirements. 

diverges at p = Pmax- For this reason, the max- 
imum density pmax cannot be exceeded. 5. The 
viscosity t] (Fig. 15) is theoretically explained. Its 
divergence at p = Pmax guarantees that velocity 
changes are smoothed out, so that the develop- 
ment of shock-like stnictures is avoided and nu- 
merical integration methods are stable. 6. The 
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Fig. 15. Representation of the theoretically ob- 
tained viscosity function. 

traffic flow is found to be stable at small and 
extreme densities, otherwise unstable (Fig. 16a). 
The result is in good agreement with Fig. 11. 
Therefore, the development of stop-and-go waves 
from almost homogeneous traffic conditions is cor- 
rectly described. The propagation velocity of the 
evolving density waves (Fig. 16b) has the right 
order of magnitude and direction, i.e. stop-and- 
go- waves move in backward direction with respect 
to the average velocity V. 



Fig. 16. The illustrations show (a) the growth rate 
A and (b) the relative backward propagation 
velocity c of small periodic disturbances with 
wave number k (corresponding to a wave 
length of i = 2-K/k), where the homogeneous 
traffic flow at density p is unstable. 

4. SUMMARY AND CONCLUSIONS 

Empirical traffic data have been compared to re- 
sults of recent macroscopic traffic models. The 
following has been found: 1. Vehicle velocities arc 
approximately normal distributed. This can be 
explained by an improved kinetic model. 2. The 
dynamics on neighboring lanes is strongly corre- 
lated so that the total freeway cross section can be 
described in an effective way. 3. The distance be- 
tween forming stop-and-go waves increases in the 
course of time, which corresponds to the merging 
of density clusters in traffic simulations. 4. The 
dynamical variance has peaks when the velocity 
changes considerably. It can be approximated by 
the variance-density relation and an additional 
term which arises from averaging the vehicle data 
over flnite time intervals. 5. The fluctuations of 
traffic density and average velocity essentially cor- 
respond to a white noise. 6. Traffic flow is imstable 
above a density of about 12 vehicles per kilometer 
and lane. This can be understood by means of 
a linear stability analysis of a new macroscopic 
traffic flow model which has been derived from the 
dynamics of driver- vehicle units via an improved 
kinetic model. 7. The model allows to reconstruct 
the velocity-density and variance-density relations 
at high densities and takes into account velocity 
fluctuations as well as the finite space require- 
ments of vehicles. 
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